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Abstract: 

We describe the development of a theoretical description of the structure of finite 
nuclei based on a relativistic quark model of the structure of the bound nucleons which 
interact through the (self-consistent) exchange of scalar and vector mesons. 

I. INTRODUCTION 

By now it is well established that one needs many-body forces to understand the 
structure of atomic nuclei. There are many ways of dealing with this problem. In the 
space available we cannot review the problem in general, rather we shall focus on re- 
cent progress based on one specific model - the quark-meson coupling model originally 
proposed by Guichon 

The quark-meson coupling model may be viewed as an extension of QHD in which 
the nucleons still interact through the exchange of a and u mesons. However, the mesons 
couple not to point-like nucleons but to confined quarks. In studies of infinite nuclear 
matter it was found that the extra degree of freedom provided by the internal structure 
of the nucleon means that one gets quite an acceptable value for the incompressibility 
once and g^^ are chosen to reproduce the correct saturation energy and density. This 
is a significant improvement on QHD at the same level of sophistication. 

In the light of current experimental work in relativistic heavy ion collisions, which 
produce nuclear matter at densities several times normal, there has been some initial 
work on the variation of baryon and meson properties with density using the quark- 
meson coupling model There have also been some interesting applications to the 
properties of finite nuclei using the local-density approximation, notably the Okamoto- 
Nolen-Schiffer anomaly |^ and super-allowed Fermi /3-decay However, the inherent 
problems of the local-density approximation mean that these applications can at best be 
semi-quantitative and it is clearly very important that the extension to finite nuclei be 
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developed. 

Our aim here is review a recently developed formulation of the quark-meson coupling 
model for finite nuclei 0, based on the Born-Oppenheimer approximation. We shall 
pay particular attention to the spin-orbit force in the model and its relation to the 
corresponding force in conventional models involving meson exchange between point-like 
nucleons. Some initial results for finite nuclei will also be presented. 

II. THE BORN-OPPENHEIMER APPROXIMATION FOR FINITE NUCLEI 

The solution of the general problem of a composite, quantum particle moving in back- 
ground scalar and vector fields that vary with position is extremely difficult. One has a 
chance to solve the particular problem of interest to us, namely light quarks confined in 
a "nucleon" which is itself bound in a finite nucleus, only because the nucleon motion is 
relatively slow and the quarks highly relativistic. Thus the Born-Oppenheimer approxi- 
mation, in which the "nucleon" internal structure has time to adjust to the local fields, is 
naturally suited to the problem. It is relatively easy to establish that the method should 
be reliable at the level of a few percent. 

Even within the Born-Oppenheimer approximation the nuclear surface gives rise to 
external fields that may vary appreciably across the finite size of the nucleon. Our 
approach has been to start with a classical "nucleon" and to allow its internal structure 
(quark wavef unctions and bag radius) to adjust to minimise the energy of three quarks 
in the ground-state of a system consisting of the bag plus constant scalar and vector 
fields, with the values at the centre of the "nucleon". (From now on we shall not put 
quotation marks on "nucleon" , but it should be remembered that our bound nucleon is a 
quasi-particle whose structure necessarily differs from that of a free nucleon.) Of course, 
the major problem with the MIT bag (as with many other relativistic models of nucleon 
structure) is that it is difficult to boost. We therefore solve the bag equations in the 
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instantaneous rest frame (IRF) of the nucleon - using a standard Lorentz transformation 
to find the energy and momentum of the classical nucleon bag in the nuclear rest frame. 

Having solved the problem using the fields at the centre of the nucleon one can 
then use perturbation theory to correct for the variation of the scalar and vector fields 
across the bag. In first order perturbation theory only the spatial components of the 
vector potential give a non-vanishing contribution. (Note that although in the nuclear 
rest frame only the time component of the vector field is non-zero, in the nucleon IRF 
there are also non- vanishing spatial components.) This extra term is a small spin orbit 
correction to the energy 

where fig is the isoscalar magnetic moment of the nucleon bag, Sg^uj is the vector potential 
felt by the nucleon with effective mass and 1]^ is a correction factor of order unity. In 
retrospect it is not surprising that the scalar magnetic moment appears, as this correction 
is associated with the effective magnetic field of the vector potential. 

The interaction in Eq.(|I]) induces a rotation of the spin as a function of time. How- 
ever, even if fig were equal to zero, the spin would rotate because of Thomas precession. 
Suppose that at time t, the spin vector is S{t) in the IRF(t). Then we expect that, at 
time t + dt the spin has the same direction if it is viewed from the frame obtained by 
boosting the IRF(t) by dv so as to get the right velocity v{t + dt). That is, the spin 
looks at rest in the frame obtained by first boosting the NRF to v{t) and then boosting 
by dv. This product of Lorentz transformation amounts to a boost to v{t + dt) times 
a rotation. So, viewed from the IRF(t + dt), the spin appears to rotate. In order that 
our Hamiltonian be correct it should contain a piece Hprec which produces this rotation 
through the Hamilton equations of motion. A detailed derivation can be found in Refs. 
and the result is 
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prec 



-V X — ■ S. 



(2) 



2' dt 

One may find the acceleration corresponding to the interaction from the Hamilton 
equations of motion. This gives 

dv 1 



-V[M*^{R) + 3glu;iR)]. 



(3) 



dt M^{R) 

If we put this result into Eq.(|D and add the result to Eq.(|l]), we get the total spin orbit 
interaction (to first order in the velocity) 



Hprec. + Hi = Vs,oXR)S ■ L, 



(4) 



where 



K.o.(^) 



2Mf{R)R _ 



^m;,{r)\ + (1 - 2f,MR)) 



(5) 



A. Centre of Mass Motion 

We have already mentioned the difficulty of boosting the bag, a problem which is 
closely related to the removal of spurious centre of mass motion. In Ref. the effective 



mass of the nucleon at each radius was computed by removing the average value of 
the square of the momentum of the three quarks, computed in the bag at each radius. 
This gives a very strong field dependence which reduces the vector potential needed to 
reproduce the correct saturation properties of nuclear matter. In Ref. we studied the 
relativistic oscillator in an external field and found that the field dependence of the cm. 
correction was, in fact, quite small. Therefore we have not followed the precription of 
Ref. [|lOl, but instead used a phenomenological cm. correction to the bag energy of the 
form —zq/Rb, which is not strongly dependent on the applied fields. As a consequence 
the vector potential in this work tends to be a little bigger (and the nucleon effective 
mass a little smaller) than in earlier work p|,|6|, p!0| , p!Tl . 



B. Quantization of the Motion of the Nucleon 



Having obtained the expressions for the energy and momentum of the bound, classical 
nucleon we can then quantize its motion. In many ways the simplest quantization proce- 
dure would be to set P — > —iWr in a non-relativistic expansion of the energy. There is 
then a small ambiguity over the ordering of V/j and M^{R) which is discussed in detail 
in Ref. [j^. An alternative procedure, which is designed to clarify the connection to QHD, 
is to quantize using the Dirac equation for the nucleon. In this case, the idea is to write 
a relativistic Lagrangian which gives equivalent expressions for the nucleon energy and 
momentum in mean-field approximation. This Lagrangian density is 

C = iip'-f ■ dip - Mn^^P + g„{a)a^'4) - g^uj^ip'^f.ip + Cmesons, (6) 

and clearly the only difference from QHD lies in the fact that the internal structure of the 
nucleon has forced a (known) dependence of the scalar meson-nucleon coupling constant, 
ga{cr) on the scalar field itself. In terms of this coupling constant the nucleon effective 
mass is 

M*^{a) = MN-g.{cr)a. (7) 

III. SELF-CONSISTENT FIELD EQUATIONS 

In the mean field approximation, the meson field operators in Eq.(|^) are replaced by 
their time independent expectation values in the ground state of the nucleus. As the 
resulting equations are in a form which closely resembles the QHD equations in Hartree 
approximation and one can relatively easily adapt existing computer programs to solve 
the quark-meson coupling model, it seems worthwhile to summarise the field equations 
here. For simplicity we shall retain only the a and uj fields, although it is straightforward 
to generalise these equations to include iso-vector mesons. 
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As explained above the nucleon satisfies the Dirac equation 



(^7-5-M^(a)-(7^7oa;)^ = 0, (8) 

where the nucleon effective mass M^{a) is given by Eq.(^) and the scalar and vector 
fields satisfy 

(-V,2 + m^)a(f) = - (J^MUa)^ {A\^m\A), (9) 
{-Wl + mlMr} = ^^(A|^V(r-)|A). (10) 



Note that the internal structure of the nucleon enters only through the scalar field 
dependence of the scalar coupling constant. In terms of the scalar charge of the nucleon 

S{r) = / duq{u — r)q{u — r), (11) 



(where q is the quark wave function in the bound nucleon), which can be expressed in 
closed form as 



we can define C{a): 



C{r) = S{r)/S{a = Q). (13) 



Then for consistency ga{(y) and C{a) must be related by 

C(a)^7.(a = 0) = -^M*(a), 



da 

It turns out that C{a) is well approximated by a linear form 



C{a) = l-ax{g,a), (15) 
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(where g^^ = (^^-(cr = 0)) so that C decreases by between 10 and 20% between free space 
and the density of normal nuclear matter 0. Indeed, in this case one can easily solve 



Eq.(|T4D for obtaining: 



N 



(16) 



In conclusion, we note for completeness the relation between the quark level coupling 
constants and those at the nucleon level 



9a = SglS{(r = 0), = Sg^. 



(17) 



IV. MORE ON THE SPIN ORBIT FORCE 

We saw earlier that the internal structure of the nucleon leads to a spin orbit coupling 
to the (isoscalar) vector potential proportional to 1 — 2/2^ (ignoring the small medium 
correction rjs). For the p meson we find the same expression but with the isovector 
nucleon magnetic moment. Now in the isoscalar case it happens that ps is approximately 
one so that 1 — 2 fig ~ — 1 which is what one obtains directly from the non-relativistic 
reduction of the Dirac equation (H). Thus one can simply use the Dirac equation without 
any serious loss of accuracy. 

On the other hand, in the isovector case one has an isovector nucleon magnetic moment 
equal to 4.7 nuclear magnetons, which is very far from unity and it appears that the 
Dirac formalism is inappropriate. However, it is well known in the one-boson-exchange 
models of the NN force, that the p coupling to the nucleon has a large anomalous piece, 
fpipcr^'^tpdj^p^. In the mean field approximation such couplings can be ignored for nuclear 
matter because the meson field is independent of position and time. The situation is 
rather different in a finite nucleus, where the time component of the vector field varies 
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with radius. In fact, in this case it is relatively straightforward to show that the non- 
relativistic reduction of the Dirac equation, including an anomalous coupling, gives a spin- 
orbit term equal to that derived in Eq.(|) provided fp/gp is chosen to be the isovector, 
anomalous magnetic moment 3.7 = (/ip — 1) — 

Clearly we could improve the accuracy of the treatment of the uj too by adding a 
small anomalous, isoscalar term with fui/duj = —0.12. It will be very interesting to 
extend these considerations to other cases - for example, the A and S hypernuclei. For 
an initial investigation of the masses of hyperons in dense nuclear matter we refer to Ref. 
1^. (Note, however, that that work used the treatment of cm. corrections to the bag 
energy which we now believe to be inappropriate - c.f. sect. 2, above.) 

V. INITIAL RESULTS 

As an initial investigation of the application of the quark-meson coupling model to 
finite nuclei we have considered the case of ^^O. For the protons one must, of course, 
include the central Coulomb repulsion. The numerical calculation was carried out using 
the techniques described by Walecka and Serot [0]. The resulting charge density for ^^O 
is shown in Fig.|l| (dotted curve) in comparison with the experimental data (hatched 
area) and QHD [|. 
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FIGURES 




FIG. 1. The charge density of O in the present model and QHD, compared with the 
experimental distribution. 

The parameters used correspond to a free bag radius of 0.8 fm - although this shrinks 
by about 2% at nuclear matter density. As the central density tended to be a little 
high in comparison with experiment we increased the model dependent slope parameter, 
a in Eq.(|l5D, by about 10% (above that calculated in the bag model) to obtain the 
results shown. The corresponding effect on the saturation energy and density of nuclear 
matter was very small. It is interesting to note, although the physical significance of the 
observation is not at all clear, that if the quark mass was taken to be around 300 MeV, 
rather than near zero (i.e. a constituent mass rather than a current quark mass) the 
density of ^^O was just as good without the need to adjust a at all. 
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VI. CONCLUDING REMARKS 



Having made so much progress in the development of the quark-meson couphng model 
there is a great deal of interest in exploring its consequences. The obvious extensions of 
the work described here and in Ref. to heavier nuclei are already underway. In view of 
the promising results for nuclear charge symmetry breaking and /5-decay obtained using 
local density approximation we are also keen to explore these applications in a genuine 
finite nucleus. For hypernuclei the natural extension (c.f. Ref. [Q) is to assume that the 
a and lo mesons couple only to the non-strange constituents. From our discussion of the 
spin orbit force in sect. 4 and the fact that the spin of the A is carried entirely by the 
strange quark, one can easily see that the A spin orbit force will arise entirely from the 
Thomas precession term. As the scalar and vector potentials tend to cancel in that term 
(c.f. Eq.(^), this means that the A spin orbit force is very naturally suppressed in this 
model - as observed experimentally. It will be important to follow this observation with 
quantitative results. 

In view of the suggestion that vector meson masses may be substantially lower in 



dense matter [|T^,|Tj] it will also be interesting to repeat our earlier work IQ with the new 
treatment of the cm. correction - i.e. with our larger scalar and vector fields. As a 
first estimate, however, we can take the lesson of Ref. that the reduction in the mass 
scales with the number of non-strange quarks and the result in the present model that 
the nucleon effective mass is of the order of 600 MeV at 2.5po to estimate that at such 
densities the effective mass of the p meson should also be around 600 MeV. This seems 
to be roughly the range needed to understand the current experiments. 

In terms of further theoretical development it will be interesting to compare the 
present model with more phenomenological, non-linear extensions of QHD - as reviewed 
recently in Ref. We would also like to consider the replacement of the effective a 
meson exchange by two-pion-exchange within a chiral quark model such as the cloudy 
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bag |jT6l . Finally one would also like to find ways to replace at least some of the repulsion 



associated with uj exchange by nucelon overlap with quark and gluon exchange. 
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